Semi-classical Theory of Conductance and Noise in Open Chaotic Cavities 
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Conductance and shot noise of an open cavity with diffusive boundary scattering are calculated 
within the Boltzmann-Langevin approach. In particular, conductance contains a non-universal 
geometric contribution, originating from the presence of open contacts. Subsequently, universal 
expressions for multi-terminal conductance and noise valid for all chaotic cavities are obtained 
classically basing on the fact that the distribution function in the cavity depends only on energy 
and using the principle of minimal correlations. 
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Transport properties of quantum systems with classi- 
cally chaotic dynamics (chaotic cavities) have recently 
become an object of intensive investigation. By far the 
most successful method to describe them is presently the 
random matrix theory (RMT) , based on assumption 
that the scattering matrix of an open chaotic cavity is a 
member of Dyson's circular ensemble of random matri- 
ces. This hypothesis was quantified to calculate aver- 
age conductance, weak localization effects, conductance 
fluctuations, shot noise, and other quantities (see |^ for 
review). It has been proven indirectly Q by showing 
equivalence with Gaussian random matrix ensembles. 

Despite this evident success, a number of questions re- 
mains unaddressed. Thus, by construction all RMT re- 
sults are universal: They only depend on the number of 
transverse channels in the leads and the time-reversal 
symmetry. It is clear, however, that certain sample- 
specific corrections exist. In particular, average conduc- 
tance must contain geometric contributions, since a cav- 
ity connected to open leads is never fully chaotic. Below 
we show that these corrections, sensitive to the shape of 
the cavity, indeed, exist, and their relative magnitude is 
of the same order as the ratio of the cross-section of the 
leads to the total surface of the cavity. 

Another problem concerns shot noise - zero-frequency 
current-current fluctuations caused by discreteness of 
electron charge [|j . For imcorrelated transfer of classical 
particles shot noise assumes the Poissonian value Sp = 
2e(/), with (/) being the average current. Fermi statistics 
induce correlations, and suppress shot noise below Sp. 
Thus, in metallic diffusive wires one has S = 2e(/)/3: 
The suppression factor F = S/Sp equals 1/3. A re- 
markable feature is that this 1/3 suppression was derived 
in two different ways, quantum- mechanically |^ (based 
on scattering approach or equivalent Green's function 
technique), and semi-classically using Boltzmann- 

Langevin approach The equivalence between these 
two is by no means evident, and the occurrence of the 
same result in the two approaches is even considered as 



a numerical coincidence |10|. Furthermore, in symmetric 



chaotic cavities (equal numbers of channels in both leads) 
RMT gives for the suppression factor |^ F = 1/4. An al- 
ternative, but still quantum-mechanical derivation of the 
same result |l^, as well as a generalization to the arbi- 
trary two-terminal case jl^,^ and multi-terminal effects 
Jl2t have been given. However, it is also desirable to de- 
velop a semi-classical analysis of shot noise in chaotic cav- 
ities, based on the Boltzmann-Langevin approach. Com- 
parison of classical and quantum results may shed new 
light onto the physics of shot noise. 

Recently, it has been realized jl^,|l^ that a model of a 
circular billiard with diffusive scattering at the boundary 
(which models surface disorder on the scale of the wave 
length) exemplifies an "extremely chaotic" system, with 
typical relaxation time of the order of the flight time, 
and, on the other hand, can be treated analytically be- 
yond RMT. Ref. [|l^ used this example to study statisti- 
cal properties of closed cavities. Below we use this model 
to study conductance and shot noise semi-classically in 
cavities connected to two open leads. In particular, we 
find corrections to the conductance which depend on the 
position of the leads, and reproduce the 1/4-suppression 
of shot noise for the fully chaotic regime. The crucial 
observation is that the distribution function of electrons 
in the cavity in the leading order is a function only of the 
energy. 

Subsequently, assuming that this is the case for a 
general chaotic cavity, and using the principle of min- 
imal correlations (representing an alternative to the 
Boltzmann-Langevin method |^), supporting the current 
conservation, we derive the expressions for the multi- 
terminal conductance and noise. It is conceptually im- 
portant that these expressions obtained semi-classically 
reproduce the RMT results available in the literature. 

Average distribution function. We consider a cir- 
cular billiard of a radius R, to which two ideal leads are 
attached, left (L) at angles 6o - (3/2 <6 <9q-\- (3/2, and 
right (R) at -a/2 < 9 < a/2 (Fig. 1). Eventually, we 
assume angular widths of the contacts, a and /?, to be 
small. In the semi-classical description, the electron is 
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characterized by the coordinate r, direction of momen- 
tum n = p/p, and energy E (which is omitted where it 
can not cause confusion). Inside the cavity the motion 
is baUistic, and the average distribution function /(r, n) 
obeys 



nVf{r,n) = 0. 



(1) 



This equation must be supplemented by boundary con- 
ditions. At the surface (denoted H.) we consider purely 
diffusive scattering for which the distribution func- 
tion of the outgoing particles is constant and fixed by 
flux conservation, 



f{r,n)^TT I {Nn') f{r,n')dn', {Nn)<0. 




FIG. 1. Geometry of a circular cavity with diffusive bound- 
ary scattering. 

Here r G f2, and N is an outward normal to the sur- 
face; we normalized dn so that J dn = 1. Furthermore, 
electrons coming from the leads are described by equi- 
librium distribution functions, fL,R{E) = 0{fiLjf — E), 
fJ-L = eV, hr = {V is the apphed voltage). We assume 
that these incoming particles are emitted uniformly to all 
directions. Thus, denoting the cross-section of the right 
lead Qr — {r = R; —a/2 < 6 < a/2}, and similarly fl^, 
we complete our system by 



f{r,n)^fL.,R, renL.R-, (iVn) < 0. 



(2) 



Due to Eq. ([^), the distribution function inside the cav- 
ity is equal to the distribution function f{9) of outgoing 
particles, taken at the appropriate point of the surface 
{r = R;6}. For this function we obtain the integral equa- 
tion 



fmn = I 
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fid'), (3) 



which has to be solved with the condition f{d)\n^n — 
Il,r [see Eq. (H)]. This equation applies to arbitrary 
sizes and positions of the leads. To progress, we assume 
the leads to be narrow (a,/3 ^ 1) and replace inte- 
grals of the type /"^^2^(6')(i6'byaF(0). We then obtain 



f(n\ "/fl + Z^A , g(0) - g(go) a/?(/3 - «) 
J yd) = m — ^ 'a 1 — , o\7 - JR) 



a + (3 ' An 
g{e)-g{0-eo) a(3 

a + 13 



Att 

with the notation 



(a + /3)2 
(/l — Jr) , 



(4) 
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The term (a//?, + /3/i?,) /{a + (3) in Eq. (|^) is uniform and 
does not depend on 9. Inside the cavity, it is responsible 
for a contribution which is position and momentum inde- 
pendent. It is also the principal term, since all others are 
proportional to the first power of a and f3. These other 
terms represent geometric corrections, which are sensitive 
to the position of the leads, and, thus, are sample-specific. 

Conductance. The average current through the cav- 
ity, / — e Jj^^ dr t;/(r, n), is expressed through the 
distribution function of outgoing particles, f{0), up to 
linear order in a and /3, 
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f{e) -Jr), (5) 



where we defined the numbers of transverse channels 
in the left and the right lead, Nl — pfRP/t''^ and 
Nr. = pFRct/irh. From Eqs. (^ and (^ we obtain the 
conductance, 
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NlNr 



2iTfi Nl + Nr 
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87r(a + (3) 
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The first term is precisely the result given by RMT, 
and the second one represents a non-universal correction, 
which depends on the position of the leads ^o- It is pos- 
itive and reaches a maximum for 9q — tt. We stress 
that this correction is not due to weak localization ef- 
fects (which can not be treated classically at all), but 
really a geometric correction, originating from the fact 
that the escape time from an open cavity is always finite, 
and the cavity is never fully chaotic. To the best of our 
knowledge, these corrections have never been discussed 

Shot noise. Now we give a semi-classical theory of 
noise suppression. We only keep the leading terms in 
a and /?, though the Boltzmann-Langevin method [^,D, 
employed below, can also give geometric corrections. We 
express the distribution function as a sum of average 
(given by Eq. (^)) and fluctuating parts. The latter one, 
Sf{x,t), obeys the Boltzmann equation, 

id/dt + VFnV)6f{x,t)=j{x,t), 

where x = {r,n, i?}. The random Langevin source j is 
zero on average, and its correlations are 



{j{x,t)j{x',t'))=iy-'S{r-r')Sit~t') 
X 6iE - E')J{r,n,n'), 



(6) 
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with V = m/{2'KlT?') being the density of states, and 

X - /') + H'„',„/'(l - /)] , 

where is the probabiUty of scattering from n 

to n' per unit time at the point r. We introduced 
the notations / = f{r,n), and /' = f{r,n'). Note 
J Jdn — J Jdn' = 0, due to the conservation of the 
number of particles. 

The probabihty W^^'{r) can be found from the fol- 
lowing considerations. It is only non-zero for nN > 0, 
n'N < 0, and under these conditions does not depend on 
n'. Thus, W^ ,^' = IWn.n where Wn,r is the probabil- 
ity of the particle n to be scattered per unit time at the 
point r. During the time interval At the particles which 
are closer to the surface then vpnNAt > are scattered 
with the probability one, and others are not scattered at 
all. Taking the limit At ^ 0, we obtain 



VFinN)S{R 
0, 



nN > 0,n'N <0 
otherwise. 



Now we define the noise spectral power at zero fre- 
quency w = in the usual way, Sij = 2/ dt{SIi{t)SIj{0)) , 
where i, j = L,R, and 6Ii(t) is the fluctuating current at 
lead i, associated with the function Sf. After standard 
transformations we obtain for S = Sub. 

S = 2e^v J dE J dndn'drTR{r,n)TH{r,n')J{r,n,n'), 

where Tji{r,n) is the probability that the particle at 
(r,n) will exit through the right lead. It obeys Eq. (|l|) 
with the boundary condition of diffusive scattering; in 
addition, it equals 1 at and at ^l- In the lead- 
ing order Tr — a/ {a -\- /3); this order, however, does not 
contribute to S, since / Jdn = / Jdn' = 0. The next 
order is that for n pointing out to the right (left) con- 
tact, Tr = 1(0). Substituting the average distribution 
function / — {a fa + (3fji)/{a + /?) into J, we obtain 
S = 2eGVF, with the noise suppression factor 



F = 



NlNr 



(7) 



coinciding with the quantum- mechanical result [TTH^ . In 
particular, for symmetric system Nl = Nr we recover 
1/4-suppression, as found previously in Ref. 

Multi-terminal conductance. The rest of the pa- 
per, motivated by Eq. (Q), we base on the assumption 
that in the leading order the distribution function in- 
side an arbitrary chaotic cavity depends only on energy, 
but not on coordinate or direction of momentum [ p^ . 
This assumption should be considered as an alternative 
to the assumption that the scattering matrix of the cavity 
is taken from Dyson's circular ensemble [^||. First we 



show that this assumption together with the requirement 
of current conservation, determines the conductance in 
the multi-terminal case. 

Consider a chaotic cavity of an arbitrary shape, con- 
nected through N ideal leads of widths Wn to the reser- 
voirs, described by equilibrium distribution functions 
,fn{E) = fpiE — eVn), fp being Fermi distribution func- 
tion. It is instructive to define the outgoing current 
Jn{E)dE through the nth lead in the energy interval be- 
tween E and E + dE, so that the total current is /„ — 
J JndE. Denoting the distribution function in the cav- 
ity as fc{E), we easily obtain Jn{E) = e^^Gn {fc - /«), 
where Gn = e^iV„/(27r?i) is the Sharvin conductance of 
the nth contact and Nn is the number of transverse chan- 
nels. Without inelastic scattering inside the cavity, the 
current in each energy interval must be conserving. From 
the condition ^ Jn = we immediately find the distri- 
bution function fc, 



/C = ^ Oinfn, Oin = Gn/ ^ G„ 



(8) 



Then, defining the conductance matrix Gmn by Im — 
Y.n GmvYn, wc obtain 



(9) 



This conductance matrix is symmetric, and for the two- 
terminal case becomes Glr = {e^ /2'kK)[NlNr/ [Nl + 
Nr)), which is the RMT result. 

Multi-terminal shot noise. The principle difficulty 
of the application of the standard Boltzmann-Langevin 
method to the noise of the chaotic cavity is that this ap- 
proach uses the collision integral explicitly. Though 
we managed above to avoid this difficulty for the circular 
cavity with surface scattering, it can hardly be overcome 
for an arbitrary ballistic cavity with chaotic dynamics. 
To resolve this problem, we formulate the principle of 
minimal correlations, which is described below. 

The actual origin of shot noise can be viewed semi- 
classically as a result of partial occupation of the states 
by electrons. In the equilibrium state with finite temper- 
ature T the partial occupation causes fluctuations of the 
distribution function with the equal time correlator 19 



{5 fix, t)Sf{x', t)) = v-H{x - x')f{x) [1 - fix)] . (10) 

The (5-function on the rhs of this equation means that the 
cross correlations are completely suppressed. We note 
that in the chaotic cavity the cross correlations should 
be also suppressed because of multiple random scatter- 
ing inside the cavity. Therefore, we assume now that 
Eq. (^|) is valid for fluctuations of the non-equilibrium 
state of the cavity, where the function fc plays the 
role of f{x) in Eq. (10). Taking into account that 
for t ^ t' the correlator obeys the kinetic equation, 
{dt + VFnV){5f{x,t)5f{x',t')) ^ |l|, we obtain the 
formula 
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x5{n-n')S{E-E')fc{l-fc), 



VFn{t - t')] 



(11) 



which describes strictly baUistic motion and is therefore 
only valid at the time scales below the time of flight. On 
the other hand, after a time of the order of the dwell time 
Td the electron becomes uniformly distributed and leaves 
the cavity through the nth contact with the probability 
Un- For times t ^ Td (which are of interest here) this 
can be described by an instantaneous fluctuation of the 
isotropic distribution 6fciE,t), which is not contained 
in Eq. The requirement of the conservation of the 

number of electrons in the cavity leads to minimal cor- 
relations between 5fc{E,t) and 5f{x,t) |^^. Thus, for 
the fluctuation of the current at the contact n we write 



5In 



epF 



dx (nNn) [Sfix, t) + SfciE, t)] , (12) 



where n„ and 7V„ denote the surface of the contact n 
and the outward normal to this contact. The require- 
ment that current is conserved at every instant of time, 
J2n^^ri = 0, eliminates fluctuations 5fc- After straight- 
forward calculations with the help of Eq. ( |lll) we arrive 
at the expression 



-2G (T + Tc), Tc 



dEfc{l~fc)- (13) 



Eq. dl^), together with Eq. (||), constitutes a finite- 
temperature multi-terminal expression for the noise 
power in a chaotic cavity. From Eq. (^ we obtain 
Tc = (e/2) a„a™(K - V^) coth (e(V; - V,n)/2T). 

At zero temperature this reproduces the noise suppres- 
sion factor (p) and the multi-terminal RMT results for 
shot noise [|l2|. At equilibrium, Smn — —^TGmn, in ac- 
cordance with fluctuation-dissipation theorem. 

In conclusion, we presented a semi-classical theory of 
transport in chaotic cavities. First, for the model of a 
circular billiard with diffusive boundary scattering, we 
reproduced the RMT results for the average conductance 
and shot noise, and, in addition, found geometrical cor- 
rections to the average conductance, which appear due to 
the fact that a cavity with open leads is not chaotic any 
more. Subsequently, we showed how the RMT results 
for multi-terminal conductance and noise in an arbitrary 
chaotic cavity may be obtained semi-classically using the 
principle of minimal correlations. 
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